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We propose to study unitary matrix ensembles denned in terms of unitary 
stochastic transition matrices associated with Markov processes on graphs. We 
argue that the spectral statistics of such an ensemble (after ensemble averaging) 
depends crucially on the spectral gap between the leading and subleading eigen- 
value of the underlying transition matrix. It is conjectured that unitary stochastic 
ensembles follow one of the three standard ensembles of random matrix theory in 
the limit of infinite matrix size TV — > oo if the spectral gap of the corresponding 
transition matrices closes slower than 1/N. The hypothesis is tested by considering 
several model systems ranging from binary graphs to uniformly and non-uniformly 
connected star graphs and diffusive networks in arbitrary dimensions. 
Submitted to Journal of Physics A; 
Version: 16th March 2001. 



I. INTRODUCTION 

The study of quantum mechanics on graphs has become an important tool for investigating 
the influence of classical dynamics on spectra, wavefunctions and transport properties of quantum 
systems. Quantum networks have been used with great success to model quantum phenomena 
observed in disordered metals and mesoscopic systems (Shapiro 1982, Chalker and Coddington 
1988); typical behaviour found in extended, diffusive systems such as localisation - delocalisation 
transitions (Freche et al 1999), multifractal properties of wavefunctions at the transition point 
(Klesse and Metzler 1995, Huckestcin and Klessc 1999), transport properties (Pascaud and Mon- 
tambaux 1999, Huckestein et al 2000) and the statistical properties of quantum spectra (Klesse and 
Metzler 1997) have been studied on graphs in the limit of infinite network size. Recently, Kottos 
and Smilansky (1997, 1999) proposed to study quantum spectra of non-diffusive graphs with only 
relatively few vertices or nodes. This work was motivated by understanding how the topology 
of the graph as well as the boundary conditions imposed at the vertices influence the statistical 
properties of the eigenvalue spectrum. A closed expression for the level spacing distribution of 
quantum graphs has recently been given by Barra and Gaspard (2000). 

The statistical properties of general quantum systems have been studied intensively over the 
last two decades or so. One finds a strong link between spectral statistics of a quantum system 
and properties of the dynamics of the underlying classical system. Numerical evidence suggests 
that the eigenvalue spectrum of a quantum system whose classical limit is chaotic follows one of 
the three standard ensembles in random matrix theory (RMT) (Mehta 1991) after suitable energy 
averaging. Quantum systems whose classical limits show mixed or integrable classical dynamics 
deviate from the RMT results. 

In this paper we will make a connection between certain unitary matrix ensembles and Markov 
processes on graphs building upon recent work by Tanner (2000) and Pakohski et al (2001). It will 
be argued that the statistical properties of an ensemble (after ensemble averaging) depend cru- 
cially on the spectrum of the transition matrix of the corresponding Markov process. Convergence 
of the ensemble average towards the three generic unitary matrix ensembles, the circular-unitary, 
-orthogonal or -symplectic ensemble (CUE, COE or CSE) is guaranteed in the limit of infinite 
matrix size only if correlations in the stochastic Markov process decay fast enough in this limit. 
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Before introducing these concepts in the following sections, we will briefly recapitulate the main 
ideas behind the quantisation of graphs. A quantum graph is essentially a network of vibrating 
strings fulfilling certain boundary conditions at the vertices. (In the usual notation of graph 
theory, we will call a directed bond between a vertex i and a vertex j of the graph an edge (ij)). 
Wave propagation on the graph is written in terms of one-dimensional plane waves moving in 
both directions along undirected edges, and one therefore considers undirected graphs in general. 
Following the notation of Kottos and Smilansky (1997), a quantum graph is characterised by an 
N dimensional, unitary matrix Ss(fc) with k being the wavenumber and N corresponds to the 
number of directed edges. It may be written in the form 

Su(*)=D(fc)V(*); (1) 

here, D is a diagonal matrix with diagonal elements da — exp(i/cLi) where £, is associated with the 
length of the edge i. The matrix elements of V contain the reflection and transmission coefficients 
for edge - edge transitions at vertices which depend on the boundary conditions chosen, see Kottos 
and Smilansky (1997, 1999) as well as Akkermans et al (2000) for details. The matrix Sb is 
a discrete quantum propagator evolving N dimensional complex wave-vectors ^ between edges 
according to f n +i = Ss(fc)^ n . The eigenfrequencies of the system are the wavenumbers k n at 
which stationary solutions exist, i.e., S(fc n ) has an eigenvalue 1. 

When studying the connection between the quantum and classical behaviour on graphs, we 
first have to define what we mean by a 'classical' dynamics on a finite network and second what 
constitutes the classical or thermodynamic limit when letting the size of the network go to infinity. 
Defining a classical deterministic dynamics on a graph which consists of a finite number of vertices 
acting as branching points of the dynamics does indeed not make sense in general (Barra and 
Gaspard 2001). Instead, one may link the probabilistic dynamics of a Markov process defined on 
the same network with the quantum evolution described by the unitary matrix S^. Kottos and 
Smilansky (1997) suggested considering the Markov process defined by the transition matrix T with 
matrix elements given by the relation = [SsiA 2 - The stochastic matrix T is itself a propagator 
describing the time evolution of a probability distribution of particles moving stochastically through 
the network with transition probabilities between adjacent edges given by the matrix elements of T. 
The degree of chaos found in such a stochastic dynamics is characterised by the decay of correlation 
of an initial probability distribution and thus by the spectrum of T. Note that the classical Markov 
process depends on the structure of the graph as well as on the boundary conditions via V but not 
on metric properties of the graph entering through the phases Lf, this is in contrast to the stochastic 
dynamics introduced by Barra and Gaspard (2001). Note also, that the transition matrices T as 
defined above describe transitions between edges of the graph, not between vertices. 

In the following we will generalise this approach by making a connection between arbitrary 
unitary matrices and their associated stochastic transition matrices. In section ||, we will briefly 
discuss the possible Markov processes which may be associated with unitary evolution on finite 
graphs and we will specify an ensemble of unitary matrices linked to a given transition matrix. We 



will in section III introduce the spectral form factor which is the statistical quantity considered 
throughout the paper, define a classical limit for networks and formulate a random matrix conjec- 
ture for unitary matrix ensembles in the limit of infinite network size. Some model systems will 



be studied in more detail in section [V 



II. UNITARY MATRICES AND ASSOCIATED TRANSITION MATRICES 



Any unitary matrix U of dimension N can naturally be linked to a stochastic transition matrix T 
of a finite Markov process with time independent transition probabilities by making the following 
connection between matrix elements of U and T, i.e. 

u ij = r ij£ * Uj — \ u ij\ — r ij ■ (2) 

The matrix T is clearly a stochastic matrix due to the unitarity of U, that is, T fulfills 
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%>0 for all i, j = 1,...N 

N 

^i y = l for alii = 1,...N. (3) 

i=i 

This implies that 1 is an eigenvalue of T, all other eigenvalues have modulus less or equal 1. The 
matrix elements t%j can be interpreted as probabilities for making a transition from a vertex i to a 
vertex j on a graph of N vertices. The structure of the corresponding graph is explicitly defined 
through T and thus U, i.e., an edge (ij) exists if and only if ty > 0. The connection (||) provides 
therefore a natural link between a unitary 'quantum' evolution and a stochastic 'classical' dynamics 
in the form 

quantum classical 
*„ +1 =U+*„ P„+i=TW (4) 

Here, ^ n is a complex N dimensional wave-vector propagating through the network and p n is 
a probability distribution; its i-th component corresponds to the probability finding a particle 
at vertex i at time n after having wandered stochastically through the network starting with a 
probability distribution po at time n = (Berman and Plemmons 1979). Note that, in contrast to 
quantum graphs briefly described in section [| the graphs defined through general unitary matrices 
may be directed. The dynamics on the graph is now defined with respect to the vertices of the 
graph, not with respect to the edges. 

It is conjectured that statistical properties of generic quantum spectra are strongly influenced 
by how fast correlations in the corresponding classical dynamics decay (Bohigas et al 1984, Berry 
1985). We may thus expect that the statistical properties of the spectrum of a unitary matrix 
U are linked to the properties of the stochastic dynamics generated by the associated transition 



matrix T. Before discussing this further in section [II, we will explore the link between unitary 
and stochastic matrices in more detail. 



A. Unitary stochastic matrices 

Not every stochastic matrix T fulfilling the conditions (^) can be associated with a unitary 
matrix as defined in (Q). Unitarity indeed requires that besides (||), also the condition 

= 1 for all j (5) 

holds. That is, T must be doubly stochastic (Marshall and Olkin 1979). One deduces immediately 
that the vector p with components pi = X/N for all i = 1, . . .N is a stationary state, i.e., it is the 
left eigenvector of T with eigenvalue 1. (It is of course also right eigenvalue of T as for all stochastic 
matrices.) This in turn implies, that the Markov process is ergodic or irreducible (Berman and 
Plemmons 1979) if the graph is connected, that is, if it is not possible to decompose the underlying 
graph into disconnected subgraphs. A general stochastic matrix is called primitive if T fe is positive 
for some k > 0, that is t^' > for all i and j (Berman and Plemmons 1979). This implies that 
the spectrum of T given in terms of the eigenvalues {A , . . . , Aat_i} with |Aj| < \Aj\ for i > j or 
in terms of the eigenexponents {A,; = — log | | } has a finite gap between the leading exponent 
Ao = and the next to leading exponents, that is, 

A = Ai - A = Ai > . (6) 

A doubly stochastic matrix is primitive if the graphs corresponding to T™ are connected for all n. 
A finite gap in the spectrum means that initial probability distributions po on the network decay 
exponentially towards the equilibrium distribution p with decay rate > Ai. 

Without going further into the theory of doubly-stochastic matrices we note that a doubly- 
stochastic matrix can be written in terms of (N — l) 2 independent parameters, say the matrix 
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elements in the first N — 1 rows and columns. These matrix elements are constrained by the 
inequalities 

N-l N-l 

^2 hj < 1 for all j; t tJ < 1 for all i; (7) 

i=i j=i 

N-l 

^2 t y - N ~ 2 > Uj >0 for all i, j 

and thus correspond to a finite domain in the (N — l) 2 dimensional parameter space. 

Not every doubly-stochastic matrix can be associated with a unitary matrix as defined in (§). The 
rows and columns of a unitary matrix have to obey orthogonality conditions which impose further 
restrictions on the matrix elements Uj. One therefore defines the subset of doubly-stochastic 
matrices T which fulfill Uj = \uij\ 2 for some unitary matrix U as unitary-stochastic transition 
matrices (Marshall and Olkin 1979). The dimension of the parameter space for unitarity-stochastic 
matrices is (JV — l) 2 as for doubly-stochastic matrices, the parameter space covered by unitary- 
stochastic matrices is, however, in general smaller than the domain specified in (^). To get precise 
bounds for the possible parameters for unitary-stochastic matrices is a non-trivial problem in 
general and is beyond the scope of this article, see Pakohski et al (2001) for details. 



B. Unitary stochastic ensembles 



Next we will focus on the space of unitary matrices related to a given unitary-stochastic matrix 
T. This space, together with a probability measure specified later, forms an ensemble which 
we will call a unitary stochastic ensemble U(N,T). These ensembles have a surprisingly simple 
structure. The number of independent parameters determining a unitary matrix uniquely is iV 2 . 
Of these (N — l) 2 parameters are fixed by the unitary-stochastic matrix T, namely the amplitudes 



ftij- After decomposing U S U(N, T) in the form 



U = DiUD 2 = 
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one finds that the remaining 2N — 1 independent parameters are the phases ipi which can take 
any values in [0, 2ir]. The building block of the ensemble is the unitary matrix U which has been 
chosen here to have real, positive matrix elements in the first row and column. It is connected 
to the transition matrix T via the relation \v,ij\ — rij. The phases of U are all contained in the 
(N — 1) dimensional complex matrix W and are fixed by the ./V(iV — 1) orthogonality conditions 
between the rows (or columns) of U. This set of equations has a discrete set of solutions for generic 
unitary stochastic matrices and N > 3 (Pakohski et al 2001), that is W and thus U is not uniquely 
determined by T. To find all possible solutions for U is a non-trivial problem in general and is 
linked to the problem of finding the maximally available parameter space for unitary stochastic 
matrices of a given dimension. 

The ensemble U(N,T) as a whole is parameterised by 2N — 1 phases <pi at most and has, for 
fixed U, the topology of a 2N — 1 dimensional torus. The symmetry properties of the ensemble 
are essentially given by the symmetries of U, we expect in particular time reversal symmetry if U 
and thus T are symmetric. The non-uniqueness will play a role only if the possible solutions U 
belong to different symmetry classes. Assuming that this is not the case, we may treat ensembles 
for fixed T but different U as equivalent and we therefore disregard the U-dependence in what 
follows. 

Taking the trivial probability measure on the parameter space </?i, . . . , ip2N-i, we can perform 
the ensemble average of a function /(U) by straightforward integration over the angles <p, that is 
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< / >U{N,T)= , , 2N _ 1 / d<Pl... dlf 2 N-lf(fl, ■ ■ ■ ,lfi2N-i;T) . (9) 

{^j Jo Jo 

The multiple-integral can be reduced to an integral over only the first N phases if / depends on 
the eigenvalues of U only. The average (^) may be written in terms of a 'time' average over an 
ergodic path on the torus. After choosing 27V — 1 rationally independent but otherwise arbitrary 
length segments Li, . . . L2N-1, one defines the trajectory 

(ipi, . . . , tp2N-i)(k) = (kLi, . . . , kL 2 N-i) mod27r 

which covers the torus uniformly when letting the fictitious time k go to infinity. The average is 
now taken over the one-dimensional parameter family 

U(fc) =Di(fc)UD 2 (fc) (10) 

where the fc-dependence in Di and D2 enters through the replacement ipi = kLi in (§)< This 
parameterisation is a generalisation of the product form in Eqn. ([!]). The average can now be 
written as 
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fco 



< / >u(n,t)= lim -r- / dkf(U(k)). (11) 

fc -»oo Ko Jo 

The importance of choosing rationally independent lengths segments L iy also stressed by Kottos 
and Smilansky (1997), becomes obvious. For rationally dependent L;'s only a lower dimensional 
subspace of the full parameter space is covered in ([ll]) which may lead to averages different from 
the full ensemble average (0). 



III. SPECTRAL STATISTICS FOR UNITARY-STOCHASTIC ENSEMBLES 

So far we proposed to divide the unitary group into unitary-stochastic ensembles (USE) which 
are defined explicitly through unitary-stochastic matrices T. We will argue now that the spectral 
statistics of unitary matrices forming a USE depend strongly on the eigenvalues of T. 

A. The spectral form factor 

In the following we identify the spectrum of a unitary matrix U of dimension N with the set of 
eigenphases {81, . . . , 9^} of U. The statistical measure used is the so called spectral from factor, 
the Fourier-transform of the spectral 2-point correlation function 

R 2 (x) = ^ < d(9)d(9 + xfd) > u(N , Tie . (12) 
d 

Here, d(9, N) = ~ @i) denotes the density of states and the mean density d is given by 

d = N/2ir (see e.g. Tanner (1999)). The average is taken over the angle 9 and a USE. After 
averaging out the 9 dependence, one recovers the Fourier coefficients in terms of the traces of U, 
i.e., one obtains for the form factor 

K(t) =< ^iTrU^I 2 >u ( n,t) (13) 

with r = n/N and the average is taken over a USE. The traces of U™ can be written as (Kottos 
and Smilansky 1997) 

(«) 

TrU" = A p eiLp ( 14 ) 
p 
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where the summation is over all periodic or closed paths of length n on the graph. Characterising 
a given periodic path by its vertex code (v\, v-i . . . v n ), Vi € {1,2,... N\ with (ujUj+i) being allowed 
transitions between vertices, one obtains, following the notation in (0), 

n n 

Lp — ^""^ (f'viVi + i i Ap — J_ J_ r ViV i+1 ■ (15) 

i=l i=l 

The form factor can thus be written as a double sum over periodic paths on the graph 

(«) 

V 



K ^ = < n E ^e^'-V) >U{N>T) (16) 



p,p' 

(«) 



5 ^TrT" + < ]T Ap^e'^-V) > a(iV)T) . (17) 



The first term in ( |l7| ) is the so-called diagonal term (Berry 1985). It stems from periodic orbit 
pairs (p,p') related through cyclic permutations of the vertex symbol code, that is, of orbits with 
vertex codes (vx, i>2 • • • v n ), (U2, V3 . . . v n , vx), ■ ■ ■ , (v n ,Vx ■ ■ ■ v n -2, v n -x)\ there are typically n orbits 
related by cyclic permutations and all these orbits have the same amplitude A and phase L. The 
corresponding periodic orbit pair contributions in (|il| ) are thus equal to which is the classical 
probability for following the given cycle for one period. Additional periodic orbit degeneracies may 
occur due to symmetries. For time reversal symmetric dynamics, for example, periodic cycles with 
symbol code (v\, V2 ■ ■ ■ v n ) and (v n , . . . v±) have identical phases and amplitudes. This leads to 
an additional symmetry degeneracy factor g in ([l7]) which is one for non-time reversal symmetric 
dynamics and two for time reversal symmetric dynamics, for example. 

The diagonal term constitutes the important connection between the form factor and the stochas- 
tic transition matrix T. The second term in (|l7| ) is a double sum over the remaining periodic orbit 
pairs. Contributions to this term which survive the ensemble average can be formulated in terms of 
periodic orbit degeneracy classes and are due to phase correlations imposed by unitarity conditions 
(Bcrkolaiko and Keating 1999, Tanner 2000). These contributions are negligible in the limit r — ► 
after ensemble averaging, but are vital to reproduce the form factor for finite r values. 

Working in the diagonal approximation valid in the asymptotic regime n — > 00 and n/N = t — » 0, 
one obtains 

K(t) wprTrT™ = gnP(ri) (18) 
where we introduce the mean return probability per vertex (Argaman et al 1993) 

P(n) = ^TrT™ . (19) 

We will argue that the spectrum of T determines whether or not the statistical behaviour of a 
unitary-stochastic ensembles U(N, T) follows RMT in the classical limit N — » 00. Before doing 
so we have to specify more precisely what we mean by the classical or thermodynamic limit of a 
stochastic dynamics on a finite graph. 



B. The classical limit and a random matrix conjecture for USE's 

In what follows we will define the classical limit of a family of stochastic Markov processes when 
letting the number of vertices and thus the dimension of T go to infinity. We thereby distinguish 
between finite systems on the one hand and extended systems on the other. A series of Markov pro- 
cesses approximating the Perron- Frobenius operator of a deterministic system acting on a bounded 
domain is a typical example of convergence to a finite classical system. The piecewise linear maps 
on the unit interval considered by Pakohski et al (2001) are particularly simple examples where 
the leading eigenvalues of the Perron- Frobenius operator are recovered already by finite transition 
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matrices. Transition matrices with increasing dimension TV resolving the phase space dynamics on 
finer and finer scales are necessary to capture more and more details of the classical dynamics for 
generic maps. We will distinguish these types of systems from extended systems consisting of net- 
works of connected, equivalent subsystems as for example the lattices shown in Fig. |^. Appropriate 
rescaling with respect to the system size is necessary here to define useful quantities describing the 
dynamical behaviour per 'unit cell'. 

To make the notion of a classical limit precise, we will adopt the following definition in what 
follows: consider a series of Markov processes given in terms of transition matrices {Ti, i = 
1, ... ,00} with Ni = dimTi < Nj = dimTj for i < j. We will say that such a series has a well 
defined classical limit corresponding to 

• a finite classical system if the integrated return probability IPi(n) — TrT™ converges uni- 
formly to a limit function IP c i(n) in the limit iVj — > 00; 

• an extended classical system if the mean return probability per vertex Pj(n) = TrT"/iVj 
converges uniformly to a limit function P c i(n) in the limit Ni — > 00. 

The semiclassical limit for a family of unitary stochastic ensembles is then defined via a family 
of unitary-stochastic transition matrices {T^} with well defined classical limit in the sense above. 

We are now able to formulate a random matrix conjecture for unitary stochastic ensembles in 
terms of the spectral gap A similar to the Bohigas-Giannoni-Schmit conjecture for general quan- 
tum systems (Bohigas et al 1984). We propose q 

The spectral statistics of a family of unitary stochastic ensembles {f/(iV,-, Tj)} with associated 
transition matrices Tj having a well defined classical limit follows one of the three random matrix 
ensembles CUE, COE or CSE in the semiclassical limit if the spectral gap A, = Aj — Aq = A^ 
ofTi decreases slower than 1/N in the classical limit, or more precisely, that there exists a Sq > 
such that 

AN- 

lim — ^ > c> for all < 5 < 5 . (20) 



The conjecture implies that USE's associated with primitive transition matrices in the classical 
limit, that is, matrices with a non-vanishing spectral gap and exponential decay of correlation, 
follow RMT-statistics. More important is, that the bound ( po|) does not exclude RMT-statistics 
for classical dynamics with algebraic decay of correlation. We also emphasis that it is the spectral 
gap which is the crucial quantity in the conjecture. No reference is made to the Kolmogorov- 
Sinai entropy or similar measures of chaos. We will pres ent systems with positive KS-entropy not 



following RMT-statistics in the classical limit in section IV. 

The 1/N - threshold condition in ( p0| ) is a consequence of the r = n/N scaling; rewriting ( 111 ) 
in the form 



\K{t)- 9 t\ 



N-l 
i=l 



A" 



< gr(N - l)e 



-XiNt 



the condition (|20|) implies that the right hand side vanishes for fixed r and N — + 00. Much less 
clear is, however, why the purely classical condition ( |20| ) implies the RMT result for K(t) for all 
r as conjectured here. This problem lies at the heart of many studies conducted in the recent past 
(Kottos and Smilansky 1999, Berkolaiko and Keating 1999, Schanz and Smilansky 1999, 2000, 



The condition that the family {T;} must have a classical limit can be relaxed. We indeed expect that 
USE's corresponding to an arbitrary series of transition matrices {Ti} with AiNi/Nf bounded from below 
by a positive constant follow RMT in the limit Ni — » 00. The limit (|2(]) is then, however, not denned in 
general and controlling the bound is difficult in practice. The series is in addition not linked to any specific 
dynamical system. 
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Tanner 2000) and will not be addressed further. We will instead consider in the next section a few 
model systems with spectral gaps both below, on and above the critical threshold and will give 
numerical results showing that the threshold condition is indeed vital for spectral statistics. 



IV. NUMERICAL RESULTS 



A. Unitary stochastic ensembles with non-vanishing spectral gap 

We will first discuss families of USE's with transition matrices having a non-vanishing spectral 
gap in the classical limit. We expect these ensembles to follow RMT-statistics for N — > oo. Two 
specific examples are considered: binary graphs with sparsely filled transition matrices and a finite 
spectral gap and fully connected graphs with uniform transition amplitudes having an infinitely 
large spectral gap. We find in both cases power law convergence of the form factor to one of the 
three RMT - ensembles. 

a. Binary graphs. We consider a special class of binary graphs with transition matrices 

_f± if j = 2imodN or j = (2i + 1) modiV , . 

13 \ 6 otherwise ' ^ ' 

where the number of vertices N is even (Tanner 2000). Every vertex has two incoming and two 
outgoing edges and the maximal number of steps to reach every vertex from every other vertex 
is log 2 N. It is easy to see that binary graphs of dimension N — 2 k , k <E IN, s o-call ed de Bruijn 



graphs (Stanley 1999), have a well defined classical limit as defined in section IIIB which is the 
dynamics of the Bernoulli shift map. All eigenvalues of the transition matrices in this family 
are zero except the leading eigenvalue A = 1; the spectral gap is infinitely large and the return 
probability IP(n) = 1 for all N = 2 k . It can furthermore be shown that every family of binary 
graphs with dimensions N = p2 k , where p > 1 is an odd integer, has a classical limit with a 
spectral gap A = log 2 independently of p. The decomposition (g|) of the USE's is unique and the 
matrices U are orthogonal consisting of N/2 nested (2 x 2) matrices of the form 

. 1 (1 1 




0.5 1 1.5 2 2.5 T 3 

FIG. 1. The form factor for binary graphs with N = 6 and N = 24; the inset shows the standard 
deviation of K(r) from the RMT result as function of N for N = 6, 8, 12, 16, 24, 32, 48, 64 and 96. 
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The spectral statistics of binary graphs in terms of periodic cycle correlations has been studied 
by Tanner (2000) and convergence towards RMT has been found considering transition matrices 
up to N = 6. The complexity of the combinatorial formulae describing the cycle length correlations 
made it difficult to advance to larger matrix sizes. The form factor can be obtained numerically by 
performing the ensemble average (|TT|) . Fig. [l] shows the ensemble averaged form factor for matrix 
sizes N = 6 and N — 24, that is, two members of the p = 3 family. Deviations from the RMT 
result for the circular unitary ensemble (CUE) clearly decrease going from N — 6 to N — 24. The 
rate of convergence is measured in terms of the mean standard deviation a averaged here over the 
t interval shown, see the inset of Fig. |l|. The numerical findings indicate a power law behaviour 
cr(N) w N~ 4 / 3 both for the p — 1 and the p = 3 family, that is, the rate of convergence does not 
depend on the size of the spectral gap. 

b. Uniformly connected star-graphs. Fully connected graphs with constant transition probabil- 
ities, that is, 



= — for all i, 7 = 1, 

1 N 



,N 



display stochastic dynamics with instant complete decay of correlations. Every initial probability 
distribution is mapped onto the equilibrium state p = (1/iV, . . . , 1/iV) in one step and all eigenvalues 
of the transition matrix are zero except for the leading eigenvalue Ao = 1. 

By interpreting the vertices i as the edges of a graph with a single central vertex, we may view 
the stochastic process as taking place on a star-shaped graph where all transitions between edges 
through the central vertex are equally likely (including the 'backscattering' processes i — > i). Note 
that binary graphs of de-Bruijn type with N = 2 k are identical to a uniformly connected star- 
graph after exactly k steps. The topological entropy h t measuring the exponential growth rate of 
periodic cycles diverges for star-graphs in the classical limit which is in contrast to binary graphs 
with h t = log 2. 

A possible choice for the matrix U in (||) defining USE's of uniformly connected transition 
matrices are symmetric Fourier matrices of the form 
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FIG. 2. The form factor for uniformly connected star-graphs and the standard deviation a. 
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A diverging topological entropy indicates a singularity in the classical dynamics and a star-graph may 
indeed serve as a model for a system with a point-like central scatterer (Berkolaiko et al 2001). 
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As a consequence of the symmetry of U, we expect COE-statistics which is indeed observed nu- 
merically, see Fig. |[ More surprisingly is the fact that convergence towards the RMT result is 
governed by the same power law as for binary graphs, see the insets of Fig. and Fig. |. One 
finds numerically a(N) « |iV -4 / 3 for star-graphs, that is, the standard deviation falls off with 
the same exponent but a slightly smaller prefactor. The rate of convergence is thus insensitive to 
the topology of the graph measured for example by the topological entropy or the spectrum of the 
transition matrix as long as a non-zero spectral gap is established. 



B. Critical ensembles and deviations from RMT-behavior 



Next, we consider two types of systems, namely quantum star-graphs and diffusive networks, for 
which the spectral gap closes exactly or faster than the critical rate A oc l/N. The deviations from 
RMT-statistics have been well studied for both types of systems but have not to our knowledge 
been considered in terms of unitary stochastic ensembles and associated transition matrices. 

c. Quantum star-graphs. Quantum star-graphs arise naturally when one quantises a graph 
with a single central vertex attached to N undirected edges. Typical boundary conditions imposed 
on the wave equation at the central vertex result in restrictions on the possible transition rates 
and backscattering is greatly favoured. The vertex scattering matrix mentioned in (|I|), which is 
essentially equivalent to the matrix U in d§), is for Neumann boundary conditions of the form 

2 

Uij = -8ij + — . (22) 

The transition matrix corresponding to this orthogonal matrix describes a Markov process of weakly 
coupled one dimensional systems with vanishing coupling strength in the limit N — > oo; the systems 
is thus extended in the sense that it consists of an increasing number of almost decoupled equivalent 
subsystems. 

The Markov processes associated with q uantum star-graphs are topologically equivalent to uni- 



formly connected graphs discussed in IV A , that is, one can move from every vertex to every other 
vertex in one step. They do, however, differ greatly in their dynamical properties; the spectrum of 
the transition matrix associated with (j|^) is highly degenerate and can be given explicitly (Kottos 
and Smilansky 1999), that is, 

4 4 

A = 0, Ai, . . . , Ajv-i = - log(l - — ) w — . 

Quantum star-graphs have therefore a critical classical spectrum with a spectral gap vanishing 
proportional to 1 /N and one finds spectral statistics intermediate between Poisson and COE / CUE 
statistics. Due to the strong enhancement of backscattering, multiple traversals of the period-1 
orbits running along the N edges give the dominant contributions to the form factor for small r = 
J* ; multiple repetitions of these orbits are invariant under cyclic permutations of the symbol code 
and the diagonal-approximation takes on a form differing from Eqn. ( |l8|) (Kottos and Smilansky 
1999), i.e., 



K(r) » P(r) = 



e 



-4t 



for fixed r and N — ► oo. The form factor approaches 1 for large r due to periodic orbit lengths 
correlations worked out in detail by Berkolaiko and Keating (1999), see also Berkolaiko et al (2001). 



d. Diffusive networks. The quantum mechanics of classical diffusive systems has been studied 
mainly in the context of Anderson localisation and insulator-metal (i.e. localisation-delocalisation) 
transitions. A variety of systems have been considered ranging from disordered conductors to 
dynamical localisation in low-dimensional Hamiltonian systems as well as various discrete network 
models, see Dittrich (1996) and Janssen (1998) for recent review articles. 
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FIG. 3. Diffusive network models in 1 and 2 dimensions where opposite sides are identified. 



The network models considered here are similar to those studied by Shapiro (1982) consisting 
of regular lattices as shown in Fig. ||. Each vertex is connected to 2d neighbouring vertices via 
undirected edges and d is identified with the dimension of the system. We choose periodic bound- 
ary conditions and consider cubical networks having the same number of vertices, L, along each 
coordinate axis. The total number of vertices is thus L d , the number of edges equals 2dL d . 

We consider a classical Markov process on these networks given by a transition matrix with 
constant transition probabilities 



U 



1 

2d 



if i and j are connected 



and L is even. The stochastic dynamics on these networks is in the continuum limit L — > oo 
in appropriately rescaled units equivalent to (i-dimensional diffusion governed by the diffusion 
equation 



(23) 



with diffusion constant D = 5-7. p(x) is the continuum limit of the discrete probability distributions 
p. The low lying eigenvalues A of Td can be recovered from solving ( p3| ) with periodic boundary 
conditions. The eigen-spectrum of (1231) is given by 



4tt 2 £>,A 
w m = JjTl^' 1 

i=l 



(24) 



where m is a (i-dimensional integer lattice vector. The eigenvalues of converge to the {uA 
spectrum in the large wavelength limit m,<L for all i. Writing the mean return probability ( pjj ) 
of the Markov process in terms of the spectrum ( f24|) , one obtains 



P(n) = ^TrTS 



m \ 



(4wDn) d /- 



J]exp 



L 2 
4Dn 



4n 2 Dn 
L 2 

d 



(25) 
(26) 



where the sums are over all lattice vectors m, 1, respectively, and the last equation is obtained by 
Poisson-summation. The return probability approaches 



P(n) = (4ttoD) 



-d/2 = ( jL\ 
\2"im J 



d/2 



for large L and n < L 2 and converges to P(n) = 1/ L d in the limit n 00 for fixed L. The system 
is thus extended having a well defined classical limit which is of course nothing but the diffusion 
process (E3I) in an infinite domain. 
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1.4 - 




FIG. 4. The form factor for one dimensional diffusive quantum networks of different size; the data are 
compared with nP(n)/2d (dashed line). 



The transition matrices corresponding to a network of vertices on a regular lattice as shown 
in Fig. ^ are not unitary stochastic. USE's on the network can be constructed when considering 
the Markov process describing transitions between adjacent edges of the network instead. Edge- 
edge transition probabilities are again chosen to be Ay = l/2d if i is an an incoming and j an 
outgoing edge meeting at the same vertex. The transition matrix of dimension N = 2 d L d is 
unitary stochastic; it is furthermore easy to see that the nonzero eigenvalues of coincide with 
the eigenvalues of T^, that is 

TrT^ = Trf n d for all n . 

The matrix consists of L coupled 2d x 2d - matrices. The corresponding unitary stochastic 
ensemble can be written, for ex ample, in terms of L coupled uniformly connected star-graphs with 



2d edges as described in [V A . Using the diagonal approximation, Eqn. (|18|), together with the 
expression for the return probability (p6|), the form factor can for small r be written as ^ 



KM « -^-7 Trf 3 = n-^P(ra) « — ^-^(Nt) 1 -^ 2 (27) 

V ; 2dL d d 2dL d y ' 2d(47rD) d / 2K ' V ; 

where we first made the diagonal approximation valid in the limit n/L d — > and n — > oo and 
the second approximation is applicable for n < L 2 . One finds in particular for one dimensional 
diffusion 



1 / N 

K(t) fts —\ — r for small r . 
K ' 2 V 2tt 

Quantum interference leads to deviations from the diagonal approximation at r w max(§^,l) 
(Schanz and Smilansky 2000) at which the form factor approaches 1, see Fig. Eigenstates of the 
USE are all localised in the classical limit, level repulsion between eigenvalues vanishes and the 
spectral statistics converges to the Poissonian limit. The form factor forms a plateau for d = 2 
and small r values, that is, 

K(t) u(n , T2) « ^ for r<^, 



3 There is yet another minor complication; the transition matrices introduced are not primitive. The 
network dynamics decomposes into two unconnected, identical networks when considering the graph for 
T|, that is, the two-step dynamics between next-to-nearest neighbour vertices and edges. The figures ^ 
and | therefore show K{t) =< ^|TrU 2n j 2 >c/(iv 1 T d ) with r = 2n/N. 
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which persists in the limit N — ► oo, see Fig. ||. We thus find a stationary distribution which 
converges neither to the CUE nor to the Poisson - limit which is a clear indication for an ensemble 
being critical. Finally, considering d > 2, we expect convergence of the form factor to the CUE 
result for large N which is confirmed by numerical calculations (not shown here). 

The influence of the dimension d on the small r behaviour of K(t) in diffusive systems has 
been described in detail by, for example, Dittrich (1996) and references there in. It is interesting 
to reconsider these results in terms of the spectral gap as described in section III B| . The spectral 



gap between the leading and next to leading eigenvalues of Td (and thus of Td) can be read off 
from Eqn. (p4|); one obtains 

_ 4tt 2 £> _ Air 2 D{2d) 2 / d 

that is, the spectral gap falls off slower than l/N for d > 2 only! The case d = 2 is critical 
and deviations from RMT indeed remain stationary in the classical limit. One-dimensional diffu- 
sion is in this sense super-critical leading to Poisson statistics. It is worthwhile noting that the 
Kolmogorov-Sinai entropy for the lattices considered is positive for all d including d — 1 and 2, that 
is, Kks = log 2d > 0. A positive KS-entropy does therefore not necessarily imply RMT-statistics. 




V. CONCLUSION 



We propose a new way to partition the unitary group into ensembles of unitary matrices specified 
in terms of unitary-stochastic transition matrices. This provides a framework to study systemat- 
ically the connections between eigenvalue statistics of unitary matrices and the properties of an 
associated classical dynamics, here the stochastic dynamics on a graph. We define a classical limit 
of a family of stochastic networks with increasing network size. This makes it possible to give a 
strict criterion distinguishing between unitary stochastic ensembles whose spectral statistics con- 
verges towards the standard RMT results for N — > oo and those whose statistics does not. Many 
questions, however, remain unanswered. The most interesting one is certainly how to link the 
universal properties of spectra of unitary matrices found on all scales (and not only for small r or 
for long-range correlations) to properties of a classical dynamics. Universality suggests a common 
principle; the attempts made so far to describe spectral properties of graphs beyond the diagonal 
approximation do, however, all rely heavily on the specific system under consideration (Kottos and 
Smilansky 1999, Berkolaiko and Keating 1999, Schanz and Smilansky 1999, 2000, Tanner 2000) 
and a general scheme is not yet in sight. One might furthermore expect that the spectral properties 
of almost all matrices within an ensemble (after applying local averaging within a given spectrum) 
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coincide with the ensemble average, an assumption which remains to be shown to be true. Finally, 
it would be interesting to study general quantum systems corresponding to a classical deterministic 
dynamics in terms of graphs. Connections between the semiclassical limit of a series of USE's and 
the semiclassical limit of a quantum map might be a way forward to understand universality in 
quantum spectra in general. 
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